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Abstract. We consider the problem of incompressible, forced, nonhelical, homo- 
O ■ geneous, isotropic MHD turbulence with no mean magnetic field. This problem is 

essentially different from the case with externally imposed uniform mean field. There 
\ is no scale-by-scale equipartition between magnetic and kinetic energies as would be 

the case for the Alfven-wave turbulence. The isotropic MHD turbulence is the end 
state of the turbulent dynamo which generates folded fields with small-scale direction 
, reversals. We propose that the statistics seen in numerical simulations of isotropic 

■ MHD turbulence could be explained as a superposition of these folded fields and 
' Alfven-like waves that propagate along the folds. 
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' The term 'MHD turbulence' embraces a number of turbulent regimes 

described by the MHD equations. Their physics can be very differ- 
ent depending on the Mach number, presence of external forcing and 
of mean magnetic field, flow helicity, relative magnitude of the ve- 
locity and magnetic-field diffusion coefficients, etc. Here we consider 
what is perhaps the oldest MHD turbulence problem dating back to 

■ Batchelor (1950): incompressible, randomly forced, nonhelical, homo- 
geneous, isotropic MHD turbulence. No mean field is imposed, so all 
magnetic fields are fluctuations generated by the turbulent dynamo. 
We are primarily interested in the case of large magnetic Prandtl num- 
ber Pr = v/r) (the ratio of fluid viscosity to magnetic diffusivity) , 
which is appropriate for the warm ISM and cluster plasmas. Numer- 
ical evidence suggests that the popular choice Pr = 1 is in many 
ways similar to the large-Pr regime. Pr 3> 1 implies that the resistive 
scale £ v ~ Pr -1 / 2 ^ is much smaller than the viscous scale t v . Thus, 
the problem has two scale ranges: the hydrodynamic (Kolmogorov) 
inertial range £q 3> t S> £ v ~ Re _3//4 ^o (£o is the forcing scale) and the 
subviscous range £ u £ ^> £•». This makes it very hard to simulate this 
regime numerically. 

For a moment, let us consider the traditional view of the fully de- 
veloped incompressible MHD turbulence in the presence of a strong 
externally imposed mean field. This view is based on the idea of Irosh- 
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Figure 1. Energy spectra in our simulations of isotropic MHD turbulence. Here 
Re A = (u 2 ) 1/2 \/v ~ Re 1/2 , where A = y/E «| Vu| 2 )/(-u 2 ))- 1/2 is the Taylor 
microscale and Re the box Reynolds number. 



nikov (1964) and Kraichnan (1965) that it is a turbulence of strongly 
interacting Alfven-wave packets. Their phenomenology, modified by 
Goldreich and Sridhar (1995) to account for the anisotropy induced by 
the mean field, predicts steady-state spectra for magnetic and kinetic 
energies that are identical in the inertial range and have Kolmogorov 
^,-5/3 sca ii n g. An essential feature of this description is that it im- 
plies scale-by-scale equipartition between magnetic and velocity fields: 
indeed, <5uk = <5Bk in an Alfven wave. Numerics appear to confirm the 
Alfvenic equipartition picture provided there is an externally imposed 
strong mean field. The reader will find further details and references in 
Chandran's review in these Proceedings. 

In the case of zero mean field, it is tempting to argue that essentially 
the same description applies, except now it is the large-scale mag- 
netic fluctuations that play the role of effective mean field along which 
smaller-scale Alfven waves can propagate. This is, indeed, what has 
widely been assumed to be true. However, the numerical simulations 
of isotropic MHD turbulence tell a very different story. No scale-by- 
scale equipartition between kinetic and magnetic energies is observed 
numerically. There is a definite and very significant excess of magnetic 
energy at small scales. This holds true both for Pr > 1 and for Pr = 1 
(figure 1). The trend is evident already at low resolutions and persists 
at the highest currently available resolution (1024 3 , see Haugen et al., 
2003 and their paper in these Proceedings). 

In order to understand what is going on, let us consider the genesis of 
the magnetic field in the isotropic MHD turbulence. As there is no mean 
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Figure 2. From kinematic dynamo to nonlinear steady state: (a) evolution of 
magnetic and kinetic energies, (b) evolution of characteristic wave numbers. Here 
k\ = ((|Vu| 2 )/(u 2 }) 1//2 = VE/X; fc rms and k\\ are defined in the text. 

field, all magnetic fields are fluctuations self-consistently generated by 
the small-scale turbulent dynamo. This type of dynamo is a fundamen- 
tal mechanism that amplifies magnetic energy in sufficiently chaotic 
3D flows with large enough magnetic Reynolds numbers (typically 
above 100) and Pr > 1 (note that both the physics and the numerical 
evidence that apply to plasmas with Pr <C 1, e.g., stellar convective 
zones and protostellar discs, are very different: see Schekochihin et al., 
2004a). The amplification is due to random stretching of the (nearly) 
frozen-in magnetic-field lines by the ambient velocity field. What sort of 
magnetic fields does the dynamo make? During the kinematic (weak- 
field) stage of the dynamo, the growth of the field is exponential in 
time (figure 2a) and the magnetic-energy spectrum is peaked at the 
resistive scale, k v ~ Pr 1 / 2 /^, and grows self-similarly (spectral profiles 
normalised by (B 2 ) are shown in figure 1; theory is due to Kazantsev, 
1968; Kulsrud and Anderson, 1992). The dynamo growth rate is of the 
order of the turnover rate of the viscous-scale eddies, which are the 
fastest in Kolmogorov turbulence. 

Thus, the bulk of the magnetic energy concentrates at the resistive 
scale. Let us inquire, however, what actually happens to the field lines? 
It turns out that the dynamo-generated fields are not at all randomly 
tangled, but rather organised in folds within which the fields remain 
straight up to the scale of the flow and reverse direction at the resistive 
scale. Figure 3b illustrates the folded structure observed numerically. 
Figure 3a explains how such fields result from repeated application 
of random shear to the field lines. The folded structure is, in fact, 
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Figure 3. The folded fields, (a) Sketch of the stretch-and-fold mechanism, (b) 
Cross-section of field strength in a simulation of nonlinear dynamo (Pr = 500). 

a measurable statistical property of the magnetic field. A variety of 
diagnostics can be used. Ott and coworkers noticed in early 1990s that 
kinematically generated fields exhibited extreme flux cancellation due 
to field reversals (see review by Ott, 1998). Our preferred description 
has been to compare the characteristic parallel and rms wave numbers 
of the field and to study the statistics of the field-line curvature (Kin- 
ney et al., 2000; Schekochihin et al., 2002a, 2004c). We have shown 
that (i) = ((|B • VB| 2 )/(S 4 )) 1 /2 ^ k u independently of Pr while 

fcrms = ((|VB| 2 ) / (B 2 )) ~ Vr l / 2 k u > fc„ (figure 2b) (ii) the bulk of the 
curvature PDF is at flow scales. This confirms the picture of direction- 
reversing fields that are straight up to the scale of the flow. Note that 
the folded structure is not static: it is a net statistical effect of constant 
restretching and partial resistive annihilation of the magnetic fields. 

One immediate implication of the folded field structure is the cri- 
terion for the onset of nonlinearity: for incompressible MHD, the back 
reaction is controlled by the Lorentz tension force B-VB ~ k\\B 2 , which 
must be comparable to other terms in the momentum equation. This 
quantity depends on the parallel gradient of the field and does not know 
about direction reversals. Balancing B-VB ~ u-Vu, we find that back 
reaction is important when magnetic energy becomes comparable to 
the energy of the viscous-scale eddies. Clearly, some form of nonlinear 
suppression of the stretching motions at the viscous scale must then 
occur. However, the eddies at larger scales are still more energetic than 
the magnetic field and continue to stretch it at their (slower) turnover 
rate. When the field energy reaches the energy of these eddies, they 
are also suppressed and it is the turn of yet larger and slower eddies 
to exercise dominant stretching. A model constructed along these lines 
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(Schekochihin et al., 2002b) leads us to expect a self-similar nonlinear- 
growth stage during which the magnetic energy grows oc t and the 
rms wavenumber of the magnetic field drops oc t -1 / 2 due to selective 
decay of the modes at the large-/c end of the magnetic-energy spectrum. 
The folded structure is preserved with folds elongating to the size £ s 
of the dominant stretching eddy. If Pr 3> \/Re, the culmination of 
this process is a situation where total magnetic and kinetic energies 
are equal, (B 2 ) ~ {u 2 ), the fields are still folded with fey ~ inverse 

box size, and k rms has dropped by a factor of Re 1//4 compared to the 
kinematic case. Since Pr 3> \/Re, we have fc rms ~ /c„Pr 1//2 Re -1 / 4 S> k v 
still below the viscous scale. Thus, the condition Pr » yTle is necessary 
for k rms and k u to be distinguishable in the nonlinear regime. It is 
possible that a second period of selective decay ensues, but this cannot 
be checked numerically because the case Pr » yTle » 1 is unresolvable 
and likely to remain so for a long time. While astrophysical large-Pr 
plasmas are certainly in this regime, the numerics can only handle 
\/Re > Pr > 1 (or Pr S> \/Re ~ 1, see Schekochihin et al. (2004c) and 
figure 3b). What can we learn from such numerical experiments? Our 
model predicts that the nonlinear-growth stage in such a situation is 
curtailed with (B 2 )/(u 2 ) ~ Pr/\/Re < 1 and magnetic energy residing 
around the viscous scale. Note that all these estimates are asymptotic 
ones, so factors of order unity can obscure comparison with simula- 
tions, in which only very moderate scale separations can be afforded. 
With this caveat, we claim that the numerically observed state where 
there is excess magnetic energy at small scales but overall (u 2 ) > (B 2 ) 
(figure 1), is consistent with our prediction. The ratio (B 2 )/(u 2 ) does 
increase with Pr, but a conclusive parameter scan is not yet possible. 
Furthermore, the reduction of fe rms in the nonlinear regime, as well as 
the elongation of the folds (decrease of k\\), are clearly confirmed by 
the numerics (figure 2b). 

Our model was based on the assumption of some effective nonlinear 
suppression of stretching motions. This does not have to mean complete 
suppression of all turbulence in the inertial range. Two kinds of motions 
that do not amplify the field can, in principle, survive. First, the velocity 
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gradients could become locally 2D and perpendicular to the field (a 
quantitative model based on such two-dimensionalisation is described 
in Schekochihin et al., 2004b and Schekochihin et al., 2004c). One 
should expect a large amount of 2D mixing of the direction-reversing 
field lines leading to very fast diffusion of the field. If this happens, 
any selective decay must be ruled out, i.e., the nonlinear-growth stage 
cannot occur and (B 2 ) cannot grow above the viscous-eddy energy. 
Since numerics support selective decay and saturated values of (B 2 ) 
are certainly far above the energy of the viscous eddies, we tentatively 
conclude that mixing due to the surviving 2D motions is not very 
efficient. The second kind of allowed motions are Alfven waves that 
propagate along the folded direction-reversing fields (figure 4). Their 

dispersion relation is = ±(bb : kk)^ 2 (B 2 } 1/2 , where b = B/B 
and k varies between the inverse length of the folds (~ box size) 
and k u (Schekochihin et al., 2002b). These waves do not stretch the 
field, do not know about the direction reversals, and have the same 
properties as standard Alfven waves. However, Fourier transforming a 
field structure sketched in figure 4 would not give Bk = Uk (scale- 
by-scale equipartition) . Instead, the magnetic-energy spectrum would 
be heavily shifted towards small scales due to the presence of direction 
reversals. The Alfven- wave component, while mixed up with folds in the 
magnetic field, should be manifest in the velocity field. We therefore 
expect the kinetic energy to have an essentially Alfvenic spectrum, 
probably fc~ 5 / 3 . This picture accounts for the excess small-scale mag- 
netic energy observed in numerical simulations. Thus, we conjecture 
that the fully developed isotropic MHD turbulence is a superposition 
of Alfven waves and folded fields. The advent of truly high-resolution 
numerical simulations should make it possible to confirm or vitiate this 
hypothesis in the near future. 

A more extended exposition of the ideas presented above, as well as 
a detailed account of the numerical evidence, are given in an upcoming 
paper by Schekochihin et al., 2004c. 
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